Back

Crop Science

Wiley

Preprints posted in the last 30 days, ranked by how well they match Crop Science's content profile, based on 18 papers previously published here. The average preprint has a 0.01% match score for this journal, so anything above that is already an above-average fit.

1
Comparison of localGEBV and Optimal Haplotype Stacking Fitness Functions using a Novel R Package: HapSelect

Shaffer, W.; Papin, V.; Carter, Z.; Brunner, S. M.; Tong, J.; Villiers, K.; Robinson, H.; Voss-Fels, K.; Hayes, B. J.; Hickey, L.; Dinglasan, E.

2026-07-13 genetics 10.64898/2026.07.08.737160 medRxiv
Top 0.1%
5.4%
Show abstract

Haplotype-based breeding strategies have emerged as promising approaches to maximize long-term genetic gain by identifying complementary parental combinations while maintaining genetic diversity. However, these methods typically require phased genotypes and more intensive workflow pipelines and skillsets. We developed a novel local genomic estimated breeding value (localGEBV) fitness function with similar intent to the optimal haplotype stacking (OHS) framework fitness function and implemented both in the novel R package, HapSelect. Our aim was to evaluate whether phased haplotypes provide additional benefit over the more easily available dosage-based unphased genotypes in highly inbred crops. A subset of bread wheat nested association mapping (NAM) population comprising 444 lines genotyped with 6,054 DArT-Seq markers was analysed. Marker effects were estimated using rrBLUP, localGEBV and haplotype effects were calculated across linkage disequilibrium-defined haploblocks, and genetic algorithms (GA) were used to identify optimal sets of 30 founders using either a localGEBV derived fitness function with unphased, dosage inputs or the OHS fitness function with phased inputs. Selected parental sets were compared with conventional truncation selection (TS) through 150 generations of forward simulation. The OHS fitness function achieved a marginally greater optimized ultimate GEBV than the localGEBV fitness function during GA optimization, with only 18 of the 30 selected founders overlapped between the two methods. Despite these differences, forward simulations demonstrated nearly identical long-term genetic gain for localGEBV and OHS-selected founders, with both approaches outperforming conventional truncation selection by maintaining greater genetic diversity and delaying the genetic plateau. The minimal difference between localGEBV and OHS is likely attributable to the high homozygosity of the population, where localGEBV and haplotype effects are nearly confounded. These results demonstrate that dosage-based localGEBV provides a practical alternative to phased haplotype approaches for parent selection in inbred crops, substantially simplifying genomic workflows while maintaining long-term breeding performance. Future work should evaluate these methods in more diverse inbred populations and outbred species, where great haplotypic diversity may increase the advantage of true haplotype-based optimizations.

2
Far-red timing uncovers cultivar-dependent yield and bolting responses in vertical-farm spinach (Spinacia oleracea L.)

McGovern, C.; Adrio, M.; Aliki, H.; Vichos, R.; Powell, W.; Sharma, R.

2026-07-13 plant biology 10.64898/2026.07.10.737849 medRxiv
Top 0.1%
5.3%
Show abstract

Far-red light (FR; 700-750 nm) is increasingly incorporated into controlled-environment lighting because it can improve photosynthetic efficiency when combined with comparatively shorter wavelengths. In long-day leafy crops such as spinach, however, FR may also promote the transition from vegetative to reproductive growth and thereby reduce marketable yield. Most studies have evaluated FR fraction, intensity or end-of-day exposure, whereas the developmental timing of FR has rarely been tested, particularly in spinach. Here, we evaluated six commercial spinach cultivars (Amador, Harp, Renegade, Responder, Rubino and Santa Cruz) in an indoor vertical farm under a common red-green-blue background (PPFD 260-264 {micro}mol m-{superscript 2} s-{superscript 1}, 12 h photoperiod, 24 {degrees}C) and four FR timing treatments: no FR (Control), FR throughout production (FullFR), FR during early development only (EarlyFR), and FR during late development only (LateFR). LateFR increased marketable fresh weight relative to Control (244 vs 224 g) and reduced flowering incidence, whereas far-red supplied during early development reduced fresh weight (158 g) and increased flowering. The magnitude of the timing response differed among cultivars: switching from EarlyFR to LateFR recovered 0 % fresh weight in Amador but 107 % in Renegade and Rubino, with the largest penalties occurring in otherwise bolt-resistant cultivars. EarlyFR also increased total chlorophyll and reduced the chlorophyll a:b ratio. These results show that FR response in spinach is strongly conditioned by developmental stage and cultivar. Although LateFR received more total far-red than EarlyFR, it behaved like the Control, indicating that the penalty was set by far-red timing rather than dose. Treatment differences in bolting and yield tracked an estimated phytochrome photostationary-state deficit during early development: a phytochrome-deficit model markedly outperformed a cumulative-dose model ({Delta}AIC = 441), and the deficit x cultivar interaction was strong (p < 0.001), with bolt-resistant cultivars losing most yield when far-red coincided with the early developmental window. We therefore propose that FR should be treated as a genotype-dependent management variable rather than as a fixed spectral input, with late application and bolt-resistant cultivars offering the most favourable combination for vertical-farm spinach production. Framed within the breeders equation, the close match between the trial and production environment and the scope for shorter breeding cycles indoors suggest that genotype and far-red timing can be optimised jointly to accelerate genetic gain.

3
Enhancing predictive accuracy of yield traits in cassava through multi-trait genomic prediction

de Freitas, G. M.; Certuche, D. S.; Jannink, J.-L.; de Oliveira, E. J.; Garcia, A. A. F.

2026-07-06 genetics 10.64898/2026.07.01.735838 medRxiv
Top 0.2%
2.3%
Show abstract

Multi-trait genomic prediction offers a practical route to improve selection for costly, complex traits in clonally propagated crops such as cassava. In a Brazilian breeding panel of 1,078 cassava clones genotyped with 25,923 SNPs and phenotyped for six agronomic traits, we compared single-trait (ST) and multi-trait (MT) GBLUP models. Stage-wise mixed models produced BLUEs that fed into ST and MT-GBLUP. We tested five cross-validation schemes that mimic breeder realities: ST baseline (CV1); naive all-traits MT prediction for unphenotyped candidates (CV2); MT prediction using auxiliary trait phenotypes in the test set (CV3); and two sparse-phenotyping regimes with missingness by trait (CV4) or by clone (CV5) at 25%, 50%, and 75% levels. The main results were that, under the ST baseline (CV1), predictive ability ranged from 0.50 for DMC and 0.45 for FRY down to 0.13 for Le.Dis. A naive full MT model (CV2) performed approximately on par with ST-GBLUP. In contrast, MT designs (CV3) that included informative auxiliary traits, such as shoot yield and combinations with plant vigor and leaf disease severity, yielded small gains for DMC with predictive ability of approximately 0.51 (+2%), while FRY predictive ability increased to approximately 0.65 (+44%), accompanied by RMSE reductions for FRY up to approximately 13.5% (e.g. RMSE approximately 6.2). Sparse-phenotyping simulations (CV4/CV5) demonstrated that MT models sustain or even improve predictive ability under realistic missing-data regimes (PA {approx} 0.62 - 0.65). Selection concordance between MT and ST top-10% sets was generally high (>0.80), and MT configurations produced measurable improvements in expected selection response and genetic gain per cycle for several target traits. These results indicate that strategically implemented MT-GBLUP, using a small set of biologically and operationally informative auxiliary traits and optimized sparse phenotyping, can materially increase predictive accuracy and selection efciency for economically critical cassava traits while reducing phenotyping burden.

4
Novel quantitative trait loci conferring broad-based resistance to root-knot nematodes in lima bean (Phaseolus lunatus)

Tajima, A. M.; Matthews, W. C.; Duong, T.; Khanh, T. D.; Baniya, A.; Penmetsa, R. V.; Parker, T.; Farmer, A.; English, S.; Diepenbrock, C.; Gepts, P.; Roberts, P. A.; Huynh, B.-L.

2026-07-09 plant biology 10.64898/2026.06.30.735594 medRxiv
Top 0.3%
1.1%
Show abstract

Lima bean (Phaseolus lunatus) is a broadly adapted, economically important leguminous crop and a susceptible host of root-knot nematodes (Meloidogyne spp.; RKN), which are a devastating plant pathogen in agricultural systems worldwide. To date, there have been few studies to elucidate the genetic determinants of RKN resistance in lima beans. Understanding the genetic mechanisms underlying resistance is essential for improving resistance traits and incorporating them into lima bean breeding programs. To assist in marker-assisted selection, we aimed to identify and map quantitative trait loci (QTLs) conferring RKN resistance-related traits. Three recombinant inbred line (RIL) populations were used in this study. Three populations were derived by crossing two RKN-resistant parents with the same RKN-susceptible parent and with each other. All populations were genotyped using genome-wide single-nucleotide polymorphism (SNP) markers. Each population was screened for root galling (RG) and RKN egg reproduction (ER) in response to M. incognita and M. javanica in greenhouse experiments. Three major QTLs were detected and mapped on chromosome Pl04 (QRk-pl04.1), Pl05 (QRk-pl05.1) and Pl10 (QRk-pl10.1) across populations. Among them, QRk-pl05.1 and QRk-pl10.1 affected levels of RG and ER of both RKN species, while QRk-pl04.1 suppressed root galling and reproduction responses of M. incognita but not of M. javanica. These chromosomal regions defined by flanking markers will help guide marker-assisted breeding and gene discovery for broad-based RKN resistance in lima beans.

5
Introducing PHJ Media: A Unique Machine Learning -Driven Basal Formulation to Overcome Recalcitrance for Multi-Genotype Micropropagation of Cannabis sativa L.

Pepe, M.; Hesami, M.; Jones, M.

2026-07-15 plant biology 10.64898/2026.07.14.738465 medRxiv
Top 0.4%
0.6%
Show abstract

Applications of tissue culture are critical for Cannabis sativa L. (cannabis), supporting clonal propagation, germplasm preservation, pathogen elimination, among other biotechnological applications. However, extensive genetic diversity associated with cannabis results in highly variable responses to in vitro conditioning, and no consensus basal media formulation exists to support reproducible micropropagation across genotypes. To address these limitations, a hybridized ensemble-NSGA-II approach was employed for concurrent optimization of individual media components to create a species specific, cultivar inclusive basal salt formulation for cannabis micropropagation. The resulting PHJ media represents a unique formulation that overcomes recalcitrance across a wide array of cannabis cultivars, facilitating improved growth and uniformity for the nine cultivars used in its development and validation. These results remain consistent from explant initiation through multiple rounds of subculture. The ability of PHJ to overcome genotypic recalcitrance is telling of its potential applicability with an array of plant species beyond cannabis. Additionally, robust performance both with and without plant growth regulators underscores the plausible use of PHJ for diverse applications beyond standard micropropagation. Ultimately, this cultivar-inclusive basal medium demonstrates utility for both scientific research and industrial-scale operations.

6
An axiomatic approach to cultivar ranking in multi-environment trials

Kondratev, A. Y.; Ianovski, E.; Voronina, E.; Crossa, J.

2026-07-01 genetics 10.64898/2026.06.27.734959 medRxiv
Top 0.4%
0.6%
Show abstract

Multi-environment trials are central to cultivar evaluation because they reveal how candidate cultivars perform across locations, years, management conditions, and stress environments. The resulting yield matrix is a rich source of data on genotype-by-environment interaction, and a wide literature on estimation, decomposition, visualisation, and prediction of yield potential and stability has flourished. However the ultimate question of which cultivar to recommend on the basis of such a matrix is often left implicit. The question is far from trivial, and in this paper we formulate cultivar recommendation as an axiomatic ranking problem. This framework is rich enough to encompass the existing literature on stability indices, as well as any other deterministic ranking procedure. We show that many commonly used stability-based procedures can violate minimal criteria of efficiency or consistency. The result of such violations is that a cultivar with uniformly high yield could be ranked below a cultivar with uniformly low yield, or the relative ranks of two cultivars could depend on whether or not a third cultivar is present in the matrix. Our results prove that under a small number of such criteria the space of admissible rules collapses to the family of power means and their limiting cases. If we further wish to allow multiplication normalisation of yield, we are left with the geometric mean as the unique solution.

7
Knowledge-guided Bayesian optimization using pre-trained LLMs speeds up the identification of superior genotypes from germplasm collection

Hamazaki, K.; Tsuda, K.

2026-07-02 bioinformatics 10.64898/2026.06.28.735149 medRxiv
Top 0.4%
0.6%
Show abstract

Background: Germplasm collections contain wide genetic diversity that is valuable for plant breeding, but conducting phenotypic evaluation for all genotypes in field trials is rarely feasible. Bayesian optimization offers a way to decide, season by season, which genotypes to cultivate in order to identify superior genotypes with fewer evaluations. However, standard Bayesian optimization commonly starts from randomly selected genotypes and mainly relies on surrogate models built from marker genotype information, while the text-based passport information that accompanies germplasm is not fully used. We examined whether pre-trained large language models can provide prior knowledge that improves these decisions in germplasm evaluation. Results: We constructed a large-language-model-guided Bayesian optimization framework that introduces large language models into two parts of the Bayesian optimization workflow. In zero-shot warmstarting, a large language model proposes initial genotypes using passport information such as cultivar name, country of origin, and subpopulation, optionally together with principal component scores derived from genome-wide single-nucleotide-polymorphism markers. In addition, we evaluated a large-language-model-based surrogate model that predicts phenotypic values for untested genotypes using in-context learning from previously evaluated genotypes. Using a rice germplasm panel and two target traits (seed number per panicle for maximization and protein content for minimization), we compared strategies. For seed number per panicle, zero-shot warmstarting with a general-purpose instruction-following model reduced the number of evaluated genotypes needed to reach the best genotype, whereas improvements were small for protein content. When genomic information was available, Gaussian-process-based Bayesian optimization was the strongest overall approach, while the large-language-model-based surrogate model outperformed random baselines and was competitive in some settings. When genomic information was not available, predictions based on passport information improved efficiency compared with fully random strategies. Conclusions: Pre-trained large language models can inject useful agronomic knowledge into Bayesian optimization for germplasm evaluation, particularly by improving early-stage genotype selection, and can also support optimization when genomic information is unavailable. As models better handle long genomic sequences together with passport information, large-language-model-guided Bayesian optimization may become a practical and explainable decision-support approach for agricultural optimization.

8
Sunrise and sunset times are the main factors that determine the flowering time of photoperiod-sensitive sorghum

Clerget, B.; Sidibe, M.; vom Brocke, K.; Raharinivo, V.; Ortiz, D.; Trouche, G.

2026-07-08 plant biology 10.64898/2026.06.12.731875 medRxiv
Top 0.5%
0.5%
Show abstract

Crop photoperiodism models assume that flowering time is primarily controlled by daylength, yet many field observations contradict this view. We previously proposed an alternative framework integrating daily changes in sunrise and sunset times (dSR and dSS). Variety trials in Madagascar and in Argentina supported this concept: mid-late sorghum varieties from the northern hemisphere flowered late or very late when sown in November and December, consistent with the higher dSR/dSS values of the southern hemisphere summer. One Malian variety, sown monthly over six years in West Africa, exhibited high interannual variability in flowering time when sown between November and February. This revealed that up to four photoperiodic responses -- two quantitative and two qualitative, occurring at different times of the year -- may coexist within a single late photoperiod sensitive variety. All responses use only dSR and dSS cues. The qualitative responses are triggered by an internal phasic coincidence, which is set by a linear relationship between dSR and dSS at the onset of plant photoperiod sensitivity, and between dSR+dSS at panicle initiation. The research model fitted data from 28 varieties grown in Mali well. It also accurately fitted the duration to PI observed in three varieties sown at tropical and temperate latitudes. HighlightThe seasonal photoperiodic adaptation of flowering time in sorghum plants may rely on several signal transduction pathways regulated by sunrise and sunset times rather than day length.

9
An in vitro regeneration system with efficient rooting in sweet orange (Citrus sinensis) supports recovery of transgenic plants

Datta, J.; Bhowmik, S. D.; Williams, B.; Kerr, S. C.

2026-07-08 plant biology 10.64898/2026.06.16.732047 medRxiv
Top 0.5%
0.5%
Show abstract

In vitro regeneration of Citrus plants is a widely used method, however, induction of adventitious roots from regenerated shoots remains a major bottleneck, limiting the recovery of healthy plants for commercial production and genomic research for crop improvement. We established an in vitro regeneration system producing profuse, healthy roots for sweet orange (Citrus sinensis cv. Benyenda) by optimising combinations and concentrations of auxins. Prior to optimising the rooting media (RTMs), we obtained a shoot regeneration rate of 90.6% from sweet orange epicotyl explants using a cytokinin, 6-benzylaminopurine (BAP). Across twelve auxin-supplemented RTMs containing different concentrations of indole-3-butyric acid (IBA) and/or 1-naphthaleneacetic acid (NAA), rooting percentages ranged from 8 - 87.5%. The combination of IBA 1.0 mg L-1 and NAA 0.1 mg L-1 promoted the best overall performance, 75 {+/-} 7.2% rooting percentage with healthy, callus-free roots ([&ge;]5 cm in length), whereas other RTMs with other auxin combinations induced callus and limited root elongation. The best-performing SRM and RTM were subsequently used for selection and recovery of transgenic sweet orange lines carrying an empty CRISPR/Cas9 construct, resulting in an 4.8% transformation efficiency. Both transgenic and non-transgenic rooted plantlets were successfully acclimatised under glasshouse conditions with a survival rate of 90%. This enhanced regeneration system overcomes rooting bottleneck and improves plant survival,enabling faster recovery of transgenic citrus lines within four months. It supports accelerated development for commercial applications and advances in citrus genetic improvement.

10
Diversity Assessment with SNP, SSR, AFLP, and RAPD Markers in Plants: A Systematic Review and Meta-Analysis

Olagunju, Y. O.; Olawuyi, O. J.

2026-07-07 plant biology 10.64898/2026.07.03.736291 medRxiv
Top 0.5%
0.5%
Show abstract

Background. DNA-based molecular markers underpin plant genetic diversity assessment, germplasm characterisation, and conservation prioritisation. Four marker systems dominate the field: Amplified Fragment Length polymorphisms (AFLPs), simple sequence repeats (SSRs), single nucleotide polymorphisms (SNPs), and random amplified polymorphic DNA (RAPDs). No quantitative meta-analysis had pooled their performance on the canonical diversity metrics: polymorphism information content (PIC), expected heterozygosity (He), and resolution power, across plants. Existing reviews are narrative, marker-restricted, or qualitatively conclusive of infeasibility. Methods. A PRISMA 2020-compliant systematic review (registered at the Open Science Framework) was executed. Eligible studies were within-study paired comparisons genotyping the same accession panel with at least two of {SNP, SSR, AFLP, RAPD} and reporting at least one diversity metric. Effect sizes were paired standardised mean differences (Hedges' g) computed under the Bernoulli-variance approximation. Random-effects REML meta-analysis used metafor 5.0.1 with Knapp-Hartung adjustment, leave-one-out, and r-sensitivity. Results. Fifteen within-study paired contrasts were eligible, distributed across three pools. Pool 2 (SSR vs SNP, He, k = 5) yielded a pooled Hedges' g of 0.494 (95% CI: -0.078 to 1.066, p = 0.075; I-squared = 90.2%; 95% PI [-0.82, 1.81]). SSRs exceeded SNPs on He in 4 of 5 studies; leave-one-out removal of the panel-size-asymmetric outlier raised the estimate to g = 0.644 (p = 0.025). Pool 3a (dominant-marker stratum, k = 6) yielded g = 0.419 (95% CI: -0.121 to 0.960, p = 0.103; I-squared = 56.5%); five of six contrasts showed SSR or AFLP exceeding RAPD on per-locus PIC. Pool 1 (PIC, k = 3, exploratory) gave a consistent direction (g = 0.453). All three pools point in the same direction: codominant or AFLP markers carry more per-locus information than the alternative being compared. Conclusions. SSR markers reported higher per-locus diversity than SNP and RAPD markers in plant within-study paired comparisons, mechanistically grounded in the SNP biallelic ceiling and the multi-allelic richness of SSRs. The effect attenuated or reversed in selfing/low-diversity panels and at the per-panel level when SNP panels exceeded approximately 1000 loci. RAPDs show the lowest per-locus information content of the four classes.

11
Multi-trait evaluation of a tomato MAGIC population identifies promising lines with improved nitrogen use efficiency (NUE)

Baraja-Fonseca, V.; Gil-Villar, D.; Bancic, J.; Renau-Morata, B.; Salud Justamante, M.; Plazas, M.; Gramazio, P.; Vilanova, S.; Perez-Perez, J. M.; Granell, A.; Molina, R. V.; Nebauer, S. G.; Prohens, J.; Arrones, A.

2026-07-15 plant biology 10.64898/2026.07.14.738388 medRxiv
Top 0.5%
0.5%
Show abstract

Nitrogen-use efficiency (NUE) is a pivotal breeding target in tomato (Solanum lycopersicum L.) to sustain production under reduced N inputs. Here, we leveraged a recently developed tomato multi-parent advanced generation inter-cross (ToMAGIC) population to identify lines with superior performance under reduced N availability. The eight founders and a core subset of 118 ToMAGIC lines were characterized with 10,684 SNP markers and evaluated under optimal (opN, 15 mM) and suboptimal (subN, 8 mM) N supply in an experiment totalling 1,576 plants, generating 48,068 data points across 61 phenotypic variables. Under both N treatments, ToMAGIC lines exhibited transgressive segregation for most traits, confirming the value of this population as a reservoir of untapped variation. Notably, under subN conditions, harvest index (Hi) increased by 29-44%, suggesting adaptive resource redistribution toward reproductive sinks. Variance partitioning revealed that agronomic and NUE-related traits were largely under genetic control, with heritability estimates frequently above 0.80 and broadly conserved across N treatments. Multivariate trait analysis identified fruit yield N concentration (NUE component, CN,y), shoot biomass N content (NAb), and shoot growth-related traits as the main drivers of treatment differentiation. Finally, proxy traits were prioritized by integrating response magnitude, heritability, trait correlations, and treatment-discriminatory power into multi-trait selection indices. This strategy generated favorable predicted genetic gains, reaching 158% for high-performance lines and 170% for subN-adapted lines, and consistently identified lines 402, 428, 518, 800, and 816 as promising pre-breeding materials. Overall, this study supports ToMAGIC as a powerful resource for developing N-efficient cultivars suited for sustainable agriculture.

12
A genetic toolkit to reduce wheat immunogenicity and incidence of celiac disease

Rottersman, M. G.; Laudencia-Chingcuanco, D.; Zhang, W.; Guzman-Lopez, M. H.; Lin, J. W.; Zhang, J.; Caseys, C.; Burguener, G.; Kim, S.; Zhang, X.; Yunusbaev, U.; Akhunov, E.; Lee, J.-Y.; Dubcovsky, J.

2026-07-08 plant biology 10.64898/2026.06.23.734071 medRxiv
Top 0.6%
0.4%
Show abstract

Celiac disease (CeD) is an immune-mediated condition triggered by wheat gluten in genetically predisposed individuals. The immune reaction in people with CeD is driven by particular gluten amino acid sequences, or immunogenic epitopes. Some of these epitopes elicit strong immune responses in the majority of CeD patients and are designated as immunodominant epitopes. Previous research has shown correlations between the amount of immunogenic wheat epitopes consumed and the onset of CeD, suggesting that reducing wheat immunogenic epitopes may reduce CeD incidence at the population level. Gluten consists of gliadins and glutenins, with gliadins having the majority of the immunodominant epitopes and glutenins playing a major role in dough strength and breadmaking quality (BMQ). This study used radiation-induced deletions, chemical mutagenesis, and natural variation in wheat (Triticum aestivum) to generate genetic stocks with reduced immunogenic epitope content. Most lines were developed in the wheat cultivar Summit, for which we produced a full genome assembly and annotation. We used exome capture to characterize these deletions and identify prolamins located within and outside the deletions. We combined different deletions and developed molecular markers to facilitate their deployment. For chromosome arms 1BS and 1DS, we generated two alternative lines: one lacking immunogenic epitopes for the development of CeD-safe genetic stocks for research purposes, and another retaining selected glutenins for breeding commercial lines with reduced immunogenicity and adequate BMQ. By making these non-transgenic genetic stocks publicly available, we aim to accelerate the development of wheat varieties with reduced immunogenicity and, eventually, a fully CeD-safe wheat.

13
Haplotypes variations of yellow stripe like (TaYSL) genes are associated with grain iron and zinc contents in wheat (Triticum aestivum L.)

Abbasi, K.; Qayyum, H.; Naseer, S.; Sun, M.; Quraishi, M. A.; Danyal, Y.; Hao, Y.; He, Z.; Rasheed, A.

2026-07-08 plant biology 10.64898/2026.06.17.732851 medRxiv
Top 0.6%
0.4%
Show abstract

The availability of pangenome and resequencing of wheat collections have facilitated the discovery of gene-trait associations in wheat. Yellow stripe-like (YSL) proteins play a key role in the uptake and translocation of metals and yet have not been fully identified and analyzed at the genome-wide level in wheat. In this study, 26 TaYSL genes were identified and divided into four distinct clades, each clade sharing similar domains and motif compositions. Most genes were upregulated under iron deficiency, whereas homoeologs of TaYSL1 were downregulated. Both SNP-based and haplotype-based association studies were used to dissect the role of TaYSLs underpinning grain iron contents (GFeC) and zinc contents (GZnC) in wheat. TaYSL6-2B and TaYSL16-1A haplotypes showed strong association with GFeC, and TaYSL14-6A showed strong association with GZnC in multiple field trials. The distribution of favorable haplotypes in global wheat collection of [~]3000 accessions showed that majority of haplotypes were more prevalent in landraces and winter wheat compared to modern cultivars and spring types, indicating their potential for use in breeding. The combination of favorable haplotypes of three YSL genes associated with GFeC and GZnC were very rare, and most of the wheat accessions has single or double favorable haplotypes. These findings provide the first comprehensive characterization of the TaYSL gene family in wheat and identify significant SNPs and elite haplotypes that can be utilized for genetic improvement and biofortification.

14
Stem photosynthesis from wild Prunus arabica enhances growth, advances bloom and increases yield in cultivated almond

Zeira, D.;Eisenbach, O.;Harel-Beja, R.;Trainin, T.;Hatib, K.;Terner, L.;Abd-Elhadi, M.;Brukental, H.;Shapira, O.;Zait, Y.;Holland, D.;Shemer, T.

2026-06-25 Plant Biology 10.64898/2026.06.23.734067 medRxiv
Top 0.6%
0.3%
Show abstract

Rising winter temperatures threaten deciduous fruit tree productivity by depleting carbohydrate reserves during dormancy. This study investigated Stem Photosynthetic Capacity (SPC), a rare adaptive trait from wild Prunus arabica, as a mechanism to enhance almond carbon economy. Using extreme segregating groups from the F1 population (P. dulcis X P. arabica), we evaluated physiological performance through high-resolution lysimetric and multi-year orchard monitoring. High-SPC [SPC(+)] genotypes maintained significantly greater stem CO2 assimilation and transpiration during leafless periods compared to low-SPC [SPC(-)] progenies. Over five successive seasons, SPC(+) trees exhibited a 33.3% increase in trunk secondary growth and reached 10% bloom approximately 8 days earlier. Most importantly, the SPC(+) group achieved a 4.6-fold increase in mean kernel yield when compared to SPC(-) group. These findings demonstrate that SPC provides a flexible, supplementary winter carbon source that directly supports both vegetative and reproductive development. Integrating SPC into commercial almond breeding programs may offer a valuable strategy to improve climate resilience and help sustain yields under warming conditions. HighlightIntegrating stem photosynthesis into commercial almond hybrids provides a winter carbon source that advances blooming, expands trunk growth by [~]33%, and increases kernel yields by more than 4.5-fold.

15
Genome-Wide Markers Predict Metribuzin Tolerance in Southern Soft Red Winter Wheat

Sellani, J.; Anzueto, H.; Arcenaux, K.; Price, P. T.; Brown-Guedira, G.; Harrison, S.; DeWitt, N.

2026-07-03 genomics 10.64898/2026.06.28.733875 medRxiv
Top 0.7%
0.3%
Show abstract

Metribuzin is a versatile herbicide effective against various annual grasses and broadleaf weeds found in wheat fields. However, it can cause foliar damage to wheat, impacting plant health and yield. A clearer understanding of the genetic architecture associated with metribuzin tolerance is necessary to guide marker-based breeding strategies. This study evaluated 351 historic Gulf Atlantic Wheat Nursery (GAWN) wheat breeding lines representative of southern US soft red winter wheat (SRWW) germplasm. Field trials were conducted at Winnsboro (WN) and Baton Rouge (BR), Louisiana, in 2016 and 2017. Metribuzin was applied at specific growth stages[DN1.1], and tolerance was assessed based on visual foliar damage. Genomic data from 6,252 filtered single nucleotide polymorphism (SNP) markers were used to estimate narrow-sense heritability, conduct genome-wide association (GWAS), and assess genomic prediction accuracy using genomic best linear unbiased prediction (GBLUP). Broad-sense heritability ranged from 0.54 to 0.69 within environments and reached 0.77 across environments, while narrow-sense heritability ranged from 0.35 to 0.47, indicating moderate additive genetic control. No SNP surpassed the significance threshold, but genomic prediction (GP) showed moderate to strong predictive ability (PA) across environments, with the highest accuracy (r = 0.62) observed between BR17 and WN17. These results indicate that metribuzin tolerance in SRWW is primarily controlled by multiple small-effect loci and that GS provides a more effective breeding strategy than marker-assisted selection for improving tolerance in southern wheat germplasm.

16
From Phenomics to Genomics: Macro-GWAS of Almond Morphology and Quality

Mas Gomez, J.; Rubio Angulo, M.; Duval, H.; Dicenta, F.; Martinez-Garcia, P. J.

2026-07-07 plant biology 10.64898/2026.07.06.736816 medRxiv
Top 0.8%
0.3%
Show abstract

In plant breeding and genetics, recent advances in high-throughput phenotyping are beginning to meet the growing demand for large-scale, high-quality phenotypic data that emerged after the development of next-generation sequencing technologies. Recent developments in phenomics have been incorporated into almond breeding programs, facilitating the large-scale acquisition of quantitative phenotypes and the dissection of the genetic architecture underlying morphological and quality-related traits. The implementation of a high-throughput phenotyping platform integrating RGB and hyperspectral imaging with genotyping using the 60K almond SNP array enabled the large-scale characterization of almond populations and the identification of 567 robust marker-trait associations across 66 traits. These analyses revealed two major genomic hotspots on chromosomes 2 and 5 associated with morphological and quality-related traits. These regions harbored biologically relevant candidate genes, including genes associated with OVATE family proteins, brassinosteroid signaling, protein ubiquitination, and acyl-CoA metabolism, as well as other regulators of organ growth, cell proliferation, hormone signaling, and seed development. Furthermore, a novel candidate gene encoding a COMT-like O-methyltransferase involved in lignin biosynthesis was identified and proposed to contribute to shell hardness, a major genetically controlled trait in almond. Together, these findings demonstrate the potential of integrating high-throughput phenomics and genomics to dissect complex traits, identify candidate genes, and accelerate genomics-informed breeding in almond.

17
SeedMeasure: an efficient approach and open-source program to quantify seed size

Sims, B.;Gaudinier, A.;Blackman, B.

2026-06-29 Plant Biology 10.64898/2026.06.27.734974 medRxiv
Top 0.8%
0.2%
Show abstract

PremiseSeed size and morphology are critical traits in agriculture, ecology, and genetics, but high-throughput quantification of these traits is often limited by labor-intensive manual measurements or expensive, platform-specific imaging software. Methods and ResultsWe developed SeedMeasure, a lightweight, open-source, and cross-platform command-line tool written in Python that automates the measurement of seed area, length, and width from images. Using a simple imaging setup, the program processes images by correcting for perspective skew, filtering debris, and exports quantitative data alongside quality-check images. We validated SeedMeasure across nine diverse species, ranging from small Arabidopsis thaliana seeds to large Zea mays kernels. The tool quickly handles images using multithreading and demonstrates high reproducibility, yielding low coefficients of variation across repeated runs. ConclusionsCompared to existing software, SeedMeasure is free, offers faster processing through parallel computing, and provides standalone executables that require no programming dependencies. SeedMeasure offers an accessible, cost-effective, and high-throughput approach for rapid phenotypic profiling, making advanced seed morphological analysis available to researchers without specialized laboratory hardware.

18
Proteomic and Metabolomic Profiling of Transgenic Pod Borer-Resistant Cowpea: Assessing Unintended Molecular Changes and Their Implications for Ecosystem Resilience

Isah, A.;Yoila, M.;Ndana, R.;Ibrahim, A.;Ogunremi, O.

2026-06-25 Plant Biology 10.64898/2026.06.24.734197 medRxiv
Top 0.9%
0.2%
Show abstract

BackgroundThe commercialization of Nigerias single-line pod borer-resistant (PBR) cowpea (IT97KT), the first transgenic cowpea variety in the world expressing Cry1Ab gene, has raised questions about potential unintended molecular changes and their ecological implications. This study employed integrated proteomic and metabolomic profiling to compare the transgenic line with its non-transgenic isoline (IT97KN) and assess molecular indicators associated with ecosystem resilience. MethodsProteomic analyses were conducted using LC-MS/MS following filter-assisted sample preparation, while metabolomic profiling employed GC-MS and UHPLC-MS/MS platforms. Differential protein and metabolite abundance were assessed using label-free quantification, volcano plot analysis, principal component analysis (PCA), hierarchical clustering, and Gene Ontology (GO) enrichment analyses. ResultsProteomic profiling revealed substantial overlap between IT97KT and IT97KN, with only a limited subset of proteins exhibiting significant differential abundance. Upregulated proteins in IT97KT were primarily associated with seed storage, redox regulation, oxidative stress mitigation, and defense-related functions, including Late Embryogenesis Abundant Protein 1 (LEA1), vicilins, thioredoxin, and iron superoxide dismutase. Among 37 proteins linked to ecological adaptation, only LEA1, CPRD22, and Bg7S showed significant differences. Similarly, only carbonic anhydrase II displayed differential abundance among proteins associated with potential ecological risk. PCA and clustering analyses demonstrated high proteomic similarity between genotypes. Metabolomic analyses identified sixteen major metabolites, predominantly fatty acids, with no statistically significant differences in abundance or composition between transgenic and non-transgenic lines ConclusionsThe transgenic PBR cowpea exhibited minimal unintended proteomic and metabolomic alterations relative to its non-transgenic isoline. These findings indicate that Cry1Ab insertion did not substantially disrupt molecular pathways associated with ecological adaptation, environmental risk, or metabolic homeostasis, providing molecular evidence supporting the environmental and biosafety equivalence of PBR cowpea.

19
Dissecting antibiosis resistance to Phthorimaea absoluta in wild and cultivated tomato accessions

Amegan, K. E.; Magot, F.; Desneux, N.; Del-Valle, S.; Salgon, S.; Kergunteuil, A.; Caromel, B.; Larbat, R.; Lavoir, A.-V.

2026-07-13 plant biology 10.64898/2026.07.11.737942 medRxiv
Top 0.9%
0.2%
Show abstract

AbstractTomato production faces a persistent challenge from the tomato leaf miner, Phthorimaea absoluta, a pest that severely limits yields while effective resistance in cultivated varieties remains scarce. To address this gap, wild tomato relatives represent a promising reservoir of resistance traits. In this study, 24 tomato accessions, including both cultivated types and wild species, were evaluated under greenhouse (no-choice) and tunnel (choice) conditions. Resistance mechanisms were characterized through measures of antibiosis such as leaflet lesion type, proportion of attacked leaflets, and mine density. The results revealed substantial variation between and within species, allowing classification of accessions into resistant, intermediate, and susceptible groups through multivariate analysis. Notably, the wild accession Solanum habrochaites PI248707 exhibited strong resistance, in contrast to susceptible cultivated varieties such as Rose de Berne. Under choice conditions, PI248707 sustained limited damage and disrupted larval development, with early instar larvae present but few reaching advanced stages, indicating an inhibitory defense response. Untargeted metabolomic profiling further highlighted pronounced constitutive differences between wild and cultivated accessions, with S. pennellii and S. habrochaites displaying higher metabolic diversity. By integrating phenotypic and metabolic data, specific metabolite classes associated with resistance were identified. These findings underscore the potential of wild tomato germplasm in breeding programs, with PI248707 standing out as a strong candidate for resistance introgression.

20
Distance-to-optimum biological drift as a new framework for interpreting routine laboratory results: a benchmark against Reference Change Values across 62 routine biomarkers

Bezier, C.; Rolland, J.; Boutin, R.; Gruson, D.

2026-07-06 biochemistry 10.64898/2026.07.06.736744 medRxiv
Top 0.9%
0.2%
Show abstract

Background: We propose the biological drift framework for the interpretation of biological test results: a z-score-like framework based on optimized and personalized reference populations and a distance-to-optimum drift metric for longitudinal interpretation relative to an estimated individual optimum. We benchmarked biological drifts against Reference Change Values (RCVs), which are used to interpret serial laboratory results by defining the minimum change expected to exceed normal within-subject biological variation CVi. Objectives: To benchmark biological drifts against the classical biological-variation framework and assess their consistency with RCV thresholds across routine biomarkers. Methods: For 62 routine biomarkers, biological drift levels were compared with RCVs after transformation to test the consistency between the two frameworks. Results: Severe biological drifts mostly exceeded the 95% RCV threshold, indicating changes unlikely to be explained by short-term biological variation alone. In contrast, moderate drifts reached the 95% RCV threshold for approximately one in two biomarkers, suggesting that many moderate distance-to-optimum deviations may remain within expected variability, particularly for biomarkers with large within-subject variation CVi. Results are particularly interesting for the follow-up of people with diabetes and for the management of thyroid and hepatic disorders. Conclusions: Biological drifts derived from optimized personalized reference populations are broadly consistent with the RCV framework for identifying biologically meaningful deviations from the optimum and may therefore be relevant for the monitoring of certain biomarkers across several medical conditions in clinical practice.